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00 ' 

f^ ' We consider reversible breaking of adhesion bonds or folding of proteins under the 

o' 

psj ■ influence of a constant external force. We discuss the stochastic properties of the 

unbinding/rebinding events and analyze their mean number and their variance in 

f^ ' the framework of simple two-state models. In the calculations, we exploit the anal- 

ogy to single molecule fluorescence and particularly between unbinding/rebinding 

^^-H , and photon emission events. Environmental fluctuation models are used to describe 

O . 

• ' deviations from Markovian behavior. The second moment of the event-number dis- 

ri ■ tribution is found to be very sensitive to possible exchange processes and can thus 

'^ . be used to identify temporal fluctuations of the transition rates. 

o 

O . PACS numbers: 82.37.Np, 82.37.Rs, 87.10.Mn, 87.15.Fh 
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The unbinding of adhesion molecules or the unfolding of proteins can be monitored with 
dynamic force spectroscopy (DFS) on a single molecule level[ll]. In many applications a 
linearly increasing load is applied to the system and the rupture force distribution is moni- 



00 

O \ tored as a function of the loading rate, yielding information about the details of the rupture 

O ' event such as the possible existence of intermediates. Similar to other single molecule tech- 



niques, such as single molecule fluorescence fSMF) or single channel recording, DFS is able 



H \ to disentangle complex reaction pathways[2|, 13|. Two setups are usually employed to de- 



a. 

termine bond-strengths or bond-lifetimes. In the linear loading mode, a force proportional 
to the observation time is applied, e.g. by moving the cantilever of an atomic force micro- 
scope away from the surface. Force-clamp spectroscopy (FCS), on the other hand, exerts 
a constant load on the bond, i.e. a fixed force is applied. In the latter case, one observes 
individual unbinding/rebinding events of the molecular system as a function of time. In 
both cases characteristic forces are in the pN-regime and bond-lifetimes are on the order 

nn 

of seconds |4, |5j. The theoretical models usually employed to describe the rupture process 
describe the dynamics in the direction of the reaction coordinate in terms of diffusive bar- 
rier crossing 6| . The rupture force or time is calculated from the mean first passage time 



in a potential consisting of the molecular potential and the harmonic potential describing 

flnri 

the application of the force [7|,|8|,|9j. So far, predominantly irreversible rupture events have 
been investigated, both theoretically and experimentally. More recently, reversible systems 



with finite rebinding rates came into focus |10l. I 111 . Il2l |. Experimentally, reversible unfold 



ing/refo! 
hairpins 



ding transitions have been observed in a variety of natural systems, such as RNA 



13 



14| and ubiquitin 151] as well as specially designed molecules such as calixarene 



catenanes[16|. In these systems fluctuations are observed due to the statistical nature of 
the transitions between a 'closed-bond' and an 'open-bond' conformation. While in the 
context of SMF it is standard to analyze the fluctuations in the photon counts a treatment 
of the mentioned fluctuations in DPS is missing. 

In the present letter we develop the theory for the analysis of the fluctuations observed 
in PCS experiments of systems that show reversible bond-breaking events. We consider 
unbinding/rebinding transitions between two states, a 'closed-bond' state A and an 'open- 
bond' state B with rates fcoff(/cxt) for A —>■ 5-transitions and konifext) for the reverse tran- 
sition, where /ext denotes the externally applied force. While for time-independent rates 
^oflf(/ext) and konifcxt), the process is Markovian, we additionally consider non-Markovian 
fluctuations. SMP has proven to be a versatile tool to observe deviations from Markovian 
behavior, for a review see [l7|. In the context of DPS a brief discussion of non-Markovian 
fluctuations has been given by Hyeon and Thirumalai[18l]. The impact of such fluctuations 
on experimental results apparently has not been considered in a quantitative manner up 
to now. However, non-Markovian fluctuations can be expected to be of importance in 
particular when considering the dynamics of protein-unfolding and refolding due to the 
complex nature of this phenomenon 19(] . 

In order to be able to describe the statistics of the unbinding/rebinding events, we 
consider the probability of finding the system in either state together with the number of 
transitions that have already occurred at time t, cf. Pig.l. 

Por the formal treatment, one considers the process consisting of the two variables 
{q{t) , m{t)} , where q{t) denotes the reaction coordinate and m{t) the number of transitions. 
In the two-state approximation the distribution is given by p{q,m;t) ~ nyi(m; t)p^(g) -t- 
nB{m;t)p^^{q) with the equilibrium populations p'xiQ), X = A, B. Here, we assumed 
that the intra- well relaxation is fast compared to the time scale of knff and kon (adiabatic 
approximation). The quantities nx{m;t) obey the master equation 20 



2q|,U: 



dtUAim; t) = -kosifext)nAifn; t) + kon{fex\)nB{m - 1; t) 



dtnB{m;t) = -kon{fext)nBi'm;t) + kos{fcxt)nAi'm;t) 



(1) 



In this form, the fluctuations in the number of transitions become evident. In particular, 
the marginal distributions nxit) = J2m=o^x{^t) and p{m;t) = Jdqp{q,m;t) can be 
computed. The global populations nx{t) are those considered in the mentioned previous 



treatments 
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U. Il2l| and nA{t) is related to the survival probability p(r) = —dnA{'T')/dT = 
fcoflf(/ext)exp(-[A;off(/ext) + fcon(/cxt)]^)- It is to be noted that g{T) = tp{t) shows strong 
deviations from Poisson behavior if /Con(/ext) ~ Ksifext), i-e. if the equilibrium constant 

^•''^-' - ti!y <^' 

is near unity. In addition, the correlation time, (r) = Jo°^^[Pi^)A ~ 

fcoflf(/ext)/ [kosifcxt) + konifcxt)] , shows a maximum as a function of the force. In the 
Bell-limit (large Ja and fe, cf. Fig.lc), k^^''^^\f) = k^^f^e'^^^ and ^n'^^H/) = ^6""^^, this 
maximum is located at f^- = In [ii'eq(0)(l + 2aB/aA)]/{c(A + C(B)- For the parameters given 
in the caption to Fig.l, one finds in this approximation f^ ~ 45.3 pN, while for the general 
Kramers rates kosif) = k^f^e^^^'^^--^^^^^ and konif) = A;°j^e"°-B/(i+///B) one has a value of 
41.4 pN, a deviation of 10%. 



The subsequent analysis of eq.([I]) proceeds in the same way as in ref.[2l|]. We introduce 
generating functions 



Gx{z;t) = J2 z"'nx{m;t) {X = A, B) 

m=0 



(3) 



from which all quantities of interest can be calculated 22] . In the present letter we focus 
on the first and second moments of the distribution of the number oi A ^^ ^-transitions, 
mt)} a.d iNim. The .at.e. will be discussed ,„ terms of Mande.s Q-pa..a.eter0. 
Q{t) = [{N{tf) - {N{t))^]/{N{t))-l. This quantity will be shown to be highly susceptible 
to subtle changes ion the energy landscape. As described in detail in ref. 21|, the generating 
functions can be used to calculate these quantities. Starting in state A initially, G{z; t) = 
Ga{z; t) + Gb{z; t) reads: 

1 



G{z;t) 



»'W+'^'«. +!£,(,;() +FW-*-=. 



2W{z) 



2W{z) 



-E2{z;t) 



(4) 



with W{z) = ^{K,^-iy + AzKe^, Ey2{z;t) = (exp [-py^i^) Jodt' kos{fext,t'))) and 
Pi/2{z) = \ [(-ft'eq + 1) T W^(-2)]- Here, we already allowed for a possible time-dependence of 
kosifext, t) due to environmental fluctuations and the brackets indicate the average over the 
disorder. However, in this case, the employed formalism can only be used if we assume a 



unique -ft'eql/ext) = kon{fcxt,t)/kos{fcxt,t) to exist. For Markovian fluctuations, one simply 
has Ei/2iz;t) = exp (-pi/2(^)A;oflf(/cxt)i)- 

The mean number of transitions {N{t)) after a short time (on the order of the inverse 
equilibrium constant), where it depends on the actual initial conditions, grows linearly 
with slope 

R = hm ^{N{t)) = M/ext)fcon(/ext) ^ ^ 

For the parameters given in the caption to Fig.l, we plot {N{t)) versus time in the up- 
per panel of Fig. 2 along with results from a kinetic Monte Carlo simulation of the same 



system 23|]. It is evident that also for times on the order of fcofr one can observe the linear 
behavior. The slope R according to eq.© is given as a function of the external force in 
the inset. It shows a maximum as a function of fext, which in the Bell-limit is located at 
fn = In [i^eq(0)(a;A/a_B)]/(«A + C(b) (deviations from Kramers results ~ 10%). 
The Q-parameter tends towards a constant value for long times: 

Qoo = lim Q{t) = -2 ^^q(^"-^) (6) 

*— ^^ [i^eq(/cxt) + if 

as shown in the lower panel of Fig.2. Again, the Monte Carlo simulation results co- 
incide with the analytical calculation, showing that Q{t) can be extracted accurately 
from noisy (experimental) data. The minimum in Qoo in the Bell-limit is located at 
fg = In [Keq{0)]/{aA + as) (deviations from Kramers results ~ 20%). 

These considerations show that it should be possible to extract both, {N(t)) and Q(t) 
from experimental data. For carrying out force-clamp experiments it is, however, highly 
advisable to carefully adjust filters that are meant to reduce instrumental noise but might 
suppress fluctuations that essentially contain valuable information about the inherent en- 
ergy landscape. The advantage of FCS is given by the fact that in addition to temper- 
ature the external force can be varied in order to adjust -ft'eq(/ext) to a prescribed value. 
Together with the information that can be gathered from the analysis of rupture force 
distributions[12], one obtains a very detailed view about the structure of the energy land- 
scape underlying the observed unbinding/rebinding transition. 

We now show that FCS also is a versatile tool to extract information about the en- 
ergy landscape in the situation where the fluctuations cannot be considered as Markovian. 
Hence, FCS can also be used to monitor the impact of environmental fluctuations or dy- 
namic heterogeneities in the reaction under consideration. To this end we consider two 
models for dynamic disorder, namely a two- configuration exchange model (TCM) and a 



model of gaussian environmental fluctuations (GM)[2l|. In the TCM it is assumed that 
for each state A and B there are two configurations [A, a), {A,b) and {B,a), {B,b). The 
exchange between configurations {X, a) and {X, b) takes place with a rate F. In a realistic 
treatment one would assume that F has a force-dependence that is different from that of the 
reaction rates. For simplicity, we neglect the force-dependency of F completely. The cor- 
rexponding reaction rates are fcoff(/cxt) for {A, a) -^ {B,a) and fcon(/ext) for the rebinding 
transition (a = a, b). We assume that k^M^t) = XaKsifext) and k^^ifcxt) = Xakonifext) 
where A;oflf(/cxt) and konifext) are the rates with the parameters given in Fig.l and Xa are 
constants. In the GM, one assumes a continuous environmental variable r{t) that under- 
goes a Ornstein-Uhlenbeck process with (force- independent) damping rate 7 and (r^) = cr^. 
One furthermore chooses kos{fext,r) = {r'^/cr'^)kofi{fext), which allows an analytical treat- 
ment, cf. ref. J24I]. In the context of SMF such gaussian models of dynamical disorder have 



been studied also including non-Markovian fluctuations 25| . The calculations of {N{t)) 
and Q{t) can be performed analytically[21] with the result that {N{t)) is not affected at all 
by environmental fluctuations in the long-time limit. For Q{t) the situation is completely 
different. For flnite F and 7, one flnds 

/-) /o (Markov) , v ^cqiJcxt) /^7^ 

where QC^^rkov) jg given by eq.([6]) and one has 

XTCM=(^'^^^^^k,Me.t)]r-' ; A'GM = (2fcoff(/ext))7~' (8) 

The important feature of eq.([7]) is that for flnite F (7) Q^o becomes positive! This is a 
clearcut indication of non-Markovian fluctuations since Qoo is always negative (or zero) for 
Markovian fluctuations. It is evident, that only for F ^ 00 (7 ^ 00), i.e. in the motional 
narrowing regime, the Markovian limit is recovered. In the limit of static disorder, F — > 
(7 -^ 0), Q{t) does not reach a limiting value in a flnite time. Instead, Q{t) increases 
linearly with time and can be written as Q{t) = Xt. The rate X is obtained from eq.(IHl) by 
replacing F (7) with unity. Moreover, eq.(IHl) shows that Q{t) does not allow to discriminate 
between different scenarios. For the two models considered only the time scale on which 
the limit is reached differs somewhat. 

In Fig. 3 we plot Qoo as a function of the external force for the GM. In order to show 
that this result also is relevant for the interpretation of experimental data, we show in 
Fig. 4 how the positive limit is reached as a function of time. Again, for not too small 7, 
Qoo is reached on the time scale of fcoff- 
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We have scrutinized the results of typical FCS experiments performed on systems show- 
ing reversible bond breakage. The fluctuations in such systems can be analyzed in much the 
same way as the photon counts in SMF. However, force spectroscopy allows to manipulate 
the reaction rates and therefore one can obtain information about the energy landscape 
that is not available from other single molecule techniques. Subtle features of the confor- 
mational dynamics displayed by reversibly bonded systems are revealed which allows to 
expand the versatility of force spectroscopy in terms of fluctuation analysis. We have also 
shown that it should be possible to discriminate between Markovian and non-Markovian 
fluctuations by examining Mandel's Q-parameter - a parameter which is experimentally 
accessible. We expect that experiments analyzed along the way discussed in this work will 
yield valuable new information regarding the details of complex chemical reactions under 
external force. 
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Figure captions 

Fig.l a): Sketch of the envisioned transitions between the bond (A) and the unbound 
state (B); 

b): A typical force-clamp trajectory from a simulation of Brownian dynamics in a 
double-well potential 

c) : Reaction rates in Kramers-approximation as function of external force; explictly, 
one has k^sif) = k^^e'^^fi^-f/^^^ and k^nif) = A;°j,e-"s/(i+///B). We choose for the 
parameters: aA/P = O.Snm, as/P = 0.7nm, /^ = SOOpN, j'b = 200pN. The rates 
for zero force are chosen as /c^g = Is^^ and /c^^ = 10*^s~^. The dotted lines are 
the Bell-limit for the rates, i.e. the limit of large j'a and Jb of the above expres- 
sions: k^^'^^^^f) = A;°ge"^-^ and k^J^^^^^f) = Ce""^^- k'^^ sets the time scale for all 
calculations. 

Fig. 2 (color online) 

upper panel: Mean number of transitions versus time for the same parameters as in 
Fig.l. The dots are results of kinetic Monte Carlo simulations. Inset: Slope R versus 
external force. 

Lower panel: Q{t) versus time for the same parameters as in Fig.l. Inset: Qoo versus 
external force. 

Fig. 3 (color online) 
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Qao versus /cxt for the same parameters as in Fig.l. The environmental fluctuations 
are modelled by the GM. 

Fig. 4 (color online) 

Q{t) versus /cxt for the same parameters as in Fig.l. The environmental fluctuations 
are modelled by the GM with 7 = 10s~^. The lines are for different -ft'cql/cxt); from 
bottom to top: 0.1, 0.2, 0.5, 1 (red), 2, 5 and 10. 
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